Magnetic field generated by r-modes in accreting quark stars 
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We show that the r-mode instability can generate strong toroidal fields in the core of accreting 
millisecond quark stars by inducing differential rotation. We follow the spin frequency evolution 
on a long time scale taking into account the magnetic damping rate in the evolution equations of 
r-modes. The maximum spin frequency of the star is only marginally smaller than in the absence 
of the magnetic field. The late-time evolution of the stars which enter the r-mode instability region 
is instead rather different if the generated magnetic fields are taken into account: they leave the 
millisecond pulsar region and they become radio pulsars. 



I. INTRODUCTION 

Since the first paper in which r-modes in rotating neu- 
tron stars were shown to be unstable with respect to the 
emission of gravitational waves [l| , it was recognized their 
relevance in Astrophysics to explain the observed distri- 
bution of the rotation frequency of stars in Low-Mass- 
X- Ray-Binaries (LMXBs) (see [2[ for a review). On the 
other hand, the instability triggered by r-modes is ac- 
tually damped to some extent by the viscosity of the 
matter composing the star: the higher the viscosity, the 
higher is the frequency of rotation of the star. This fact 
opens the possibility of investigating the internal compo- 
sition of compact stars by studying the viscosity of the 
different possible high density phases which can appear 
in these stellar objects. A number of papers are cur- 
rently present in the literature about the shear and bulk 
viscosities of nucleonic matter [3|- 7| , hyperonic matter [Sj- 
Il4l . kaon condensed matte r lla . Il6j , pure quark phases 
[17h27| |. and mixed phases [281]. Interestingly, by study- 
ing the so called "window of instability" of the r-modes, 
which is determined by the shear and the bulk viscosity 
of the matter, it was pointed out in [295 that a detection 
of a sub-millisecond rotating star would indicate the exis- 
tence of a very viscous phase in the star with only quark 
or hybrid stars as possible candidates. 

Also the temporal evolution of the spin frequency of a 
star under the effect of r-modes instability has been stud- 
ied for different types of composition and different phys- 
ical systems: newly born compact stars and old stars in 
binaries [29l - [3^ . While neutron stars could be very pow- 
erful gravitational waves emitter only during their first 
months after birth [sij . hyperonic and quark or hybrid 
stars could turn into steady gravitational waves sources 
if present in LMXBs [32, 33, 35]. 

R-modes are also responsible for differential rotation in 
the star which in turn generates a toroidal magnetic field: 
besides the viscosity of matter, the production of this 
magnetic field represents a very efficient damping mech- 
anism for the r-modes. This effect has been proposed 
and investigated in 36n3S] for the case of a newly born 
neutron star and it has been recently included within the 
r-mode equations of the neutron stars in LMXB [H, : 



by calculating the back-reaction of the magnetic field on 
the r-modes instability, it has been proven that magnetic 
fields of the order of 10^^ G can be produced. Remark- 
ably, this mechanism could be at the origin of the enor- 
mous magnetic field of magnetars. 

In this paper we extend the calculations of [s^ H^l 
to the case of quark stars in LMXB. The motivation 
for this investigation is that the evolution of accreting 
stars and their internal magnetic field depends strongly 
on the r-modes instability window which, for quark stars, 
is qualitatively different from the one of neutron star. 
Indeed, due to the large contribution to the bulk viscos- 
ity of the non leptonic weak decays occurring in strange 
quark matter, the window of instability splits into two 
windows, a small one at large temperatures (which is ac- 
tually irrelevant for the evolution) and a big one at low 
temperatures [ll,|2§]. Moreover, quark stars do not have 
a crust (or just a very thin crust) and, as we will discuss, 
this prevents the trapping of the internal magnetic field 
developed during the evolution with possible observable 
signatures. 

The paper is organized as follows: in Sec. II we intro- 
duce the system of equations which provide the temporal 
evolution of a compact star in a LMXB. In Sec. Ill we 
show the results of our numerical calculations and finally 
in Sec. IV we present of conclusions. 



II. THE TEMPORAL EVOLUTION MODEL 

A. R-modes equations 

Let us review the derivation of the equations for the 
evolution of a compact star in presence of the r-niodcs. 
Following the discussion of Ref. [33], we start from the 
conservation of the angular momentum. The total an- 
gular momentum of the star can be decomposed into an 
equilibrium angular momentum J* and a canonical angu- 
lar momentum Jc^ proportional to the r-modes amplitude 
a: 
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where J» = 7*0, with fl the angular velocity and 7* = 
IME? the moment of inertia of a star of mass M and 
radius R. For a n=l polytrope / = 0.261, iiTc = 9.4 x 
10~^ (see Ref. 'Sl']) and TCj is of the order of unity (the 
results of our calculations are rather insensitive to the 
value of Kj). 

The equations for the evolution of a star are based on 
two simple considerations: 

• the canonical angular momentum, proportional to 
the r-modes amplitude, increases with the emission 
of gravitational waves and decreases due the pres- 
ence of damping mechanisms (as e.g. viscosity and 
the torque due to an internal magnetic field B): 
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where TGw{M,rt) and Tdamp{M,Q,T, B) are the 
gravitational wave emission timescale and the 
damping timescale, respectively, where the latter 
includes all the mechanisms muffling the r-modes. 
Finally with T we mean the average temperature 
of the star. 

• the total angular momentum takes contributions 
from the mass accretion and decreases due to the 
emission of gravitational waves and electromag- 
netic waves. The electromagnetic wave emission is 
associated to the presence of an external poloidal 
magnetic field which is not aligned with the rota- 
tion axis. Indicating with the braking timescale 
due to the external magnetic field, the second equa- 
tion reads: 
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where J a is the variation of the angular momentum 
due to the mass accretion. Following Ref. [s^, we 
assume Ja = M{GMRY/^. 

Combining together the equations ([T]), ([2]) and it 
is easy to obtain the evolution equations for i7 and for 
the r-modes amplitude a: 
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Here we adopt the estimate given in Ref. Q for the 
gravitational radiation rate due to the I = m = 2 current 
multipole: 
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where we have used the notation Mi. 4 
= i?/10 km and P_3 = P/l ms. 
Following Ref. [40] , the total damping rate is given by 
the sum of the shear viscosity and bulk viscosity rates 
plus the damping rate due to the internal magnetic field: 

^damp = + + ^m] (7) 

The estimates of the viscosity dampin g ra tes for the case 
of pure quark matter are given in Ref. [32| . For the shear 
viscosity it is found: 



Ts = 3.4 X W'al/'M-j/'R\l^'T^/' s 
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where as is the strong coupling constant and Tg = T/IQ^ 
K. For the bulk viscosity, the viscosity coefficient takes 
the form 



(9) 



where uj is the frequency of the r-modes in the corotating 
frame and a and /3 are coefficients given in Ref. (l7j . Due 
to this behavior, bulk viscosity is very large at T ~ 10® if 
and gets weaker at both lower and higher temperatures. 
The instability window is then splitted into a Low Tem- 
perature Window (LTW) (at T < 10® K) and an High 
Temperature Window (HTW) (T > 10® K). Due to the 
fast cooling rate of quark stars, the (HTW) is crossed 
very rapidly and does not have a big impact on the evo- 
lution of the star. So we are interested only in the LTW. 
The bulk viscosity damping timescale in the low temper- 
ature limit (T < 10® K) is given hy^: 

tI"^ ^7.9Ml,R^o'P-3Ti'm^o% ^ (W) 

where mioo is the mass of the strange quark in units of 
100 MeV. Notice that the value bulk viscosity damping 
timescale is strongly dependent on the value of mioo. 

The expression of the magnetic damping rate has 
been derived in [s^, H^, where it has been shown 
that while the star remains in the instability region, 
the r-modes generate a differential rotation which can 
greatly amplify a pre-existing magnetic field. More 
specifically, if a poloidal magnetic field was originally 
present, a strong toroidal field is generated inside the 
star. The energy of the modes is therefore transferred to 
the magnetic field and the instability is damped. 
We assume that the internal stellar magnetic field B is 
initially based on the solution obtained by Ferraro [4l| 
which can be written 14011 
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where Bd is the strength of the equatorial magnetic field 
at the stellar surface. 

To estimate the magnetic field produced by r-modes we 
start by writing the I ^ m = 2 contribution to the per- 
turbation velocity: 
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Following Ref. [37[ we get the total azimuthal displace- 
ment from the onset of the oscillation at io up to time i, 
which reads: 

Ai'^(r,i)= [ 5v'l'{t')dt' 

J to 

^ \ (^) ^'^^^ /* o?{m{t')dt' + Oia') (13) 

where k2{e) = (l/2)^(5!/7r)(sin26l - 2cos26i). The re- 
lation between the new and the original magnetic field 
inside the star in the Lagrangian approach reads [STj : 
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This equation implies that the radial dependence of the 
initial and final magnetic field is the same. Integrating 
on time the induction equation in the Eulerian approach 
one gets [ssj : 
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where B^ is the toroidal component. 

The expression of the magnetic damping rate reads pol |: 
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where E is the energy of the mode, Em is the magnetic 
energy and A « 0.99. The time integral over the r-mode 
amplitude a takes contribution from the period during 
which the star is inside the instability region. 



B. Temperature evolution 

The equations derived above are strongly dependent on 
the value of the temperature, therefore we need to com- 
pute also the thermal evolution of the star. For the sake 
of simplicity we assume that the temperature is uniform 
in the star and we make use of the estimates derived in 



Ref. [3^, for the case of quark stars. The heat capacity 
reads: 



Cv = 1.5 X lO^^M^^^RioTs erg/K 
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The star cools down essentially by emitting neutrinos, 
through the so called URCA processes. Taking into ac- 
count the direct URCA process, which are likely to occur 
in the core of quark stars, the cooling rate due to emission 
of neutrinos reads: 13211: 
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If the star is accreting mass then one has to consider the 
heating due to the accreted material. In particular, for a 
quark star the heating source comes from the conversion 
of nucleons into strange matter. Estimating the energy 
release from this conversion to be around 20 MeV per 
nucleon [s^ li^l , the heating rate due to the accretion is 
given by: 
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where Af_8 is the accretion rate in units of lO^^M© per 
year. The last contribution to the heating of the star 
comes from the viscosity. Taking into account the con- 
tributions of both bulk and shear viscosity, the associated 
heating rate reads |3 1| : 

Ey^scos^tv = 2a^nHlR^j{- + -) crg/s (20) 

where the expression for Tg and are given in eqs. ([5]) 
and (HOl), respectively, and J = (2/3) x 9.4 x lO^^/. 

Finally, the equation for the thermal evolution of the 
star is given by: 

— EdQQ^^^liQYi ~\- EyigQQgily En,eutri7io{^^^ 

In the absence of r-modes, the temperature is deter- 
mined by a balance between accretion heating and neu- 
trino cooling: 



T ^ 8.25 X lO^M^sRwY'^K 



III. RESULTS AND DISCUSSION 



(22) 



By solving the Eas. (|4l5ll6l21l) it is possible to obtain 
the temporal evolution of quark stars, taking into ac- 
count the formation of toroidal magnetic fields due to 
r-mode instability. In particular, in the present work we 
show the results only for the case of quark stars inside of 
LMXBs. Eq. (IT5|) gives the strength of the new toroidal 
component created by the winding up of a pre-existent 
poloidal field, due to the differential rotation induced by 
r-modes. Several issues remain open concerning the way 
the new magnetic fields, produced by the damping of the 
r-modes, are affected by possible instabilities. In the sta- 
bly stratified matter of a stellar interior there are two 
types of instabilities: the Parker (or magnetic buoyancy) 
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FIG. 1. The quark star instability window for three values 
of the strange quark mass. The dashed lines indicate the 
evolution of a star in a LMXB during the accretion stage up 
to the moment in which the star enters the instability window. 



and the Tayler instabilities (or pinch- type), both driven 
by the magnetic field energy in the toroidal field. The 
condition for the Tayler instability to set in is given by 

mill: 

where cja = B / (AnpY^'^r is the Alfven frequency, A^^ ~ 
5 X 10*s~^ is the compositional contribution to the buoy- 
ancy frequency and ~ 10~ ^ f {ms)T^cm^ s^^ is the mag- 
netic diffusivity which can be obtained from the electri- 
cal conductivity a by using the relation r] — l/(/iocr) 
[45I and f{ms) takes into account the dependence of 
the electron fraction on the mass of the strange quark 
and it ranges from 11 to 4 to 2.5 when changes from 
100 to 200 to 300 MeV. From Eq. ^ we can conclude 
that in a pure quark star, the Tayler instability sets in 
for B > Bcr 10^^ G. After the development of the 
Tayler instability, the toroidal component of the field 
produces, as a result of its decay, a new poloidal com- 
ponent which can then be wound up itself, closing the 
dynamo loop. When the differential rotation stops the 
field can evolve into a stable configuration of a mixed 
poloidal-toroidal twisted-torus shape inside the star with 
an approximately dipolar field connected to it outside the 
star [4g - |4g |. Concerning the buoyancy instability, it oc- 
curs for Bcr ^ 10^^ G [50], therefore it is not relevant in 
our calculations. 

It is very important to point out that the formation and 
the further stabilization of magnetic fields in the core of 
compact stars, leads to different scenarios depending of 
the star composition. For instance, neutron stars contain 
a highly conductive crust acting as a screen for the inter- 
nal large magnetic field. Conversely, whether quark stars 
contain or not a crust is still controversial. Various types 
of crust have been suggested: a tenuous crust made of 
ions suspended on the electric field associated with the 
most external layer of a quark star [l^l ; a crust made of 
a mixed phase of electrons and of quark nuggets [5l|; a 



FIG. 2. Instability window and trajectories of the evolution of 
the frequency and temperature of the star for different values 
of the mass accretion rate. The short thick lines represent 
the paths followed by the star after the instability windows is 
touched. The dots indicate the onset of the Tayler instability 
at which the evolution is stopped. 



crust made of a mixed phase of hadronic and quark mat- 
ter [sl] . The electric properties of the crust area strongly 
connected with the electron density, which is suppressed 
in the core (what has been taken into account above, 
while discussing the magnetic diffusivity) and which can 
be enhanced in the crust, at least in the model in which a 
fraction of hadron is present in the crust [52]. Although 
detailed calculations of the structure of the crust of a 
quark star are not yet available, we can assume that if 
a crust exists its electric properties interpolate between 
those of a normal neutron star and those of a bare quark 
star, the two limiting situations. 

As discussed above, when Tayler instability sets in, a 
poloidal component with strength similar to the toroidal 
one {B ~ 10^^ G) is generated. If the crust is not 
present (the case with a crust will be shortly discussed 
later in the conclusions) such large poloidal component is 
not screened and diffuses outside of the star, preventing 
the star from further accreting mass. Therefore, follow- 
ing this scenario, if a quark star inside of a LMXB pro- 
duces, due to the r-mode instability, an internal magnetic 
field larger than the Tayler instability threshold, then it 
should stop accreting mass. Thus this represents, at least 
in principle, an additional constraint on the highest ro- 
tational frequency of the compact stars in LMXBs. 
There are many parameters that enter into the equations 
for the temporal evolution of the frequency of the star. 
We fix the mass and the radius of the star to be IAMq 
and 10 km respectively. A crucial quantity for the bulk 
viscosity is the mass of the strange quark for which we 
will consider the three values nis = 100 — 200 — 300 MeV. 
The smallest value for is basically the current mass 
of the strange quark as used within MIT bag-like mod- 
els and the last value is compatible with results of chiral 
quark models as the NJL model. The astrophysical input 
parameters concern the mass accretion rate, for which we 
use the values in the range (10~^" — 10~*) MQyr~^, and 
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FIG. 3. Temporal evolution of a for different values of the 
mass accretion rate. Here the initial poloidal magnetic field 
is Bd = 10** G and = 100 MeV. 



the initial dipolar magnetic field varied in the interval 
Bo = 10^ - 10^ G. 

Let us start our discussion with the r-modes instability 
window shown in Fig. [TJ With increasing values of mg, 
the bulk viscosity increases and therefore the instability 
region becomes smaller. Also shown by dashed lines are 
the trajectories of the evolution of a star in a LMXB 
during the accretion stage and until the star enters the 
instability window for two extreme values of the mass ac- 
cretion (10~^° and 10"® M^yr'^ ). 
The evolution is calculated by starting from a configu- 
ration below the instability window: an initial value of 
100 — 200 Hz is taken as initial frequency, a is set to zero 
and the initial temperature is the equilibrium tempera- 
ture as expressed by Eq. ip^ . Due to accretion, the star 
is spun up until it reaches the instability window and at 
that point the r-modes instability starts to develop. The 
bulk viscosity dissipates part of the r-mode energy into 
heat with a consequent reheating of the star and the in- 
ternal magnetic field starts to grow. At the same time, 
the frequency of the star continues to grow due to the 
accretion torque, and the amplitude of the r-modes in- 
creases. The simultaneous effect of reheating and accre- 
tion leads the star to follow the border of the instability 
window. This is clearly shown in Fig. [51 where the thick 
lines indicate the paths followed by the star (for different 
values of the accretion rate) and the dots, at which the 
evolution stops, signal the onset of the Tayler instabil- 
ity. The temporal evolution of the r-modes amplitude is 
shown in Fig. [3] and Fig.|4l where we take rag = 100 MeV 
and TTT-s = 200 MeV, respectively. For both the cases the 
initial dipolar magnetic field has a value of 10* G. 

Notice that after an initial stage in which a oscillates 
(the star enters and exits the instability region), a then 
increases steadily but its value is still so small that r- 
modes cannot reach the saturation regime. So our cal- 
culations are independent of the value of a in the sat- 
uration regime. Parallel to the evolution of a we show 
in Figs. [5] and [5] the temporal evolution of the internal 
magnetic field. As expected, the value of follows the 
same evolution of a. Again, we stopped the evolution as 
soon as reaches the threshold for the Tayler instability 
Bcr 10^2 G. 



FIG. 4. As in Fig. |3]with = 200 MeV. 



le+12 
le+10 
le+08 
le+06 
10000 
100 



m^=100 MeV 
B =10"g 



— lO"" M^^ 


/year 


- io^'m^^ 


/year 


- 10-'° M 


^^year 



1000 10000 
time[yrsj 



FIG. 5. Temporal evolution of the internal toroidal magnetic 
field for different values of the mass accretion rate. Here 
Bd = 10* G and m, = 100 MeV. 



As discussed in Ref. [321, which the effect of inter- 
nal magnetic damping was not considered, a quark star 
inside a LMXB can be spun up up to a maximum fre- 
quency, corresponding to the frequency at which the ac- 
cretion torque balances the spin down torque due to the 
emission of gravitational waves. However, taking into ac- 
count the internal magnetic damping, a quark star with- 
out a crust should stop accreting after the onset of Tayler 
instability is reached. Therefore an important question 
is how much the limiting frequency is reduced with re- 
spect to the case in which the internal magnetic field is 
not taken into account. To this purpose, we computed 
the maximum frequencies for both the cases, estimat- 
ing in this way the effect of the internal magnetic field. 
These frequencies are shown in Fig. [7] as a function of rus 
and for two values of the initial poloidal magnetic field. 
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FIG. 6. As in Fig. [S]with rUs = 200 MeV. 
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FIG. 7. Maximum spin frequencies of quark stars as a func- 
tion of the mass of the strange quark. The results obtained 
assuming a initial poloidal magnetic field Ba = 10* G (dashed 
line) and Bd = 10® G (dot-dashed line) are compared to the 
result obtained without the internal magnetic field (solid line) . 
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FIG. 8. Trajectories in the plane P-P of quark stars after 
the development of the Tayler instability. The green line cor- 
responds to a quark star without a crust; the blue line is 
obtained if in the crust the main diffusive process is the am- 
bipolar diffusion. 



Bd = 10« G (dashed line) and = 10^ G (dot-dashed 
line) . Notice that the effect of the internal magnetic field 
is rather small, reducing the maximum rotational fre- 
quency only by a few tens Hertz. Moreover, the proven 
existence of stars in LMXBs rotating at frequencies larger 
than 600 Hz rules out a value of 100 MeV. 

Finally, we investigate the possible evolutionary scenar- 
ios of quark stars beyond the onset of the Tayler instabil- 
ity. Let us first consider a star without a crust: as soon 
as the Tayler instability sets in, the new magnetic con- 
figuration prevents the star from further accreting mass. 
The new poloidal component, of the same order of mag- 
nitude of the toroidal component ('^ lO^'^ G), will act 
as a strong braking torque, and the star will lose angu- 
lar momentum. Such evolutionary path is plotted in the 
plane P-P in Fig.([S]) and is depicted by the green line. 
It is very interesting to notice that, starting from the re- 
gion of LMXBs, the star evolves into the region of radio 
pulsars. 

On the contrary, a highly conductive crust could screen 
the internal magnetic field for a very long time. A possi- 
ble example of such a crust is that described in Ref. [11] 
where the most external layer is made of an admix- 
ture of hadrons and quarks. In that case the electron 
fraction (and therefore also the electric conductivity) is 
large in the crust and decreases towards to core of the 
star. The evolution of the star depends on the dom- 
inant dissipation mechanism associated with the crust. 
The most relevant one is probably the ambipolar diffu- 
sion (see Refs. ^,1131), whose typical timescale is: 

taruMp ~ 3 X lO^^^yr (24) 

, where L5 = L/10^ cm is the size of the region embed- 
ding the magnetic field. Due to this diffusion mechanism. 



the internal magnetic field can "quite rapidly" (in a few 
millions years) diffuse outside of the crust, and also in 
this case the star should stop accreting. In the mean- 
while, the external magnetic field should spin down the 
star. In Fig. [5] we indicate with a blue line a sketch of the 
path followed by the star when the ambipolar diffusion is 
present; also in this case the star evolves into the region 
of the radio pulsars. 

It is important to remark that the time spent by 
the quark star after the Tayler instability and before it 
reaches the region of the radio pulsars is of the order of 
a few million years, to be compared with the time spent 
inside the radio pulsar region, which is of the order of 
a few hundreds million years. Therefore we can expect 
that only a few percent of the stars following the trajec- 
tory indicated in Fig([51) will be detected before reaching 
the radio pulsar region. 

IV. CONCLUSIONS 

We have shown that strong toroidal fields can be gen- 
erated in the core of an accreting millisecond quark star 
which enters the r-mode instability window. Tayler in- 
stability sets in when the generated toroidal fields exceed 
the critical value B^^^ ~ 10^^ G and a new poloidal com- 
ponent of similar strength is produced. Our results show 
that the maximum spin frequency for quark stars does 
not change significantly when taking into account the in- 
ternal generated magnetic fields. 

The scenario after the development of the Tayler insta- 
bility depends on the presence and the properties of a 
possible crust. If the crust is not present, the gener- 
ated large poloidal component diffuses quickly outside 
the core and prevents the further accretion of mass on 
the star. On the other hand, if a highly conductive crust 
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is present, it could screen to some extent the internal 
magnetic field. However, taking into account the am- 
bipolar diffusion, which is, in this case, the dominant 
dissipation mechanism, the star could expel the internal 
magnetic field in a few millions years, which would then 
stop the accretion. In both cases the quark star evolves 
into the region of radio pulsars, as shown in FiglS] this 
represents a new possible scenario for the formation of 
radio pulsars. 



ACKNOWLEDGMENTS 



L.B. is supported by CompStar a research program of 
the European Science Foundation. G.P. is supported by 
the German Research Foundation (DFG) under Grant 
No. PA1780/2-1. J. S. B. is supported by DFG through 
the Heidelberg Graduate School of Fundamental Physics. 



N. Andersso n, [Astrophys. J.| 502, 708 (1998), 

[^Xl^:gr-qc/97 06075] 

N. Andersson an d K. D. Kokkotas, |Int. J. Mod. Phys.| 

DIO, 381 (2001), arXiv:gr-qc/0010102' 

R. F. Sawyer, Phys. Rev. D39, 3804 (1989) 

P. Haensel and R. Scliaeffer, P hys. Rev.| D45, 4708 

(1992). 

P. Haensel, K. P. Levenfish, and D. G. Yakovlev, Astron. 
Astrophys. 357, 1157 (2000), |arXiv:astro-p h/0004183 
P. Haensel, K. P. Levenfish, and D. G. 
Yakovlev, Astron. Astr ophys. | 327, 130 (2001), 

arXiv:astro- ph/010 3290 

O. Benhar and M. Valh, Phys. Rev. Lett. 99, 232501 



(2007), a rXiv:0707.2681 [nucl-th]| 

L. Lindbl om and B. J. Owen, |Phys. Rev] D65, 063006 



(2002), arXiv:astro- ph/0110558r 
P. Haensel, K. P. Levenfish, and D. G. 
Yakovlev, Astron. Astr ophys.] 381, 1080 (2002), 

arXiv:astro-ph/ 0110575 

D. Chatterjee a nd D. Ba ndyopadhyay, jPhys. Rev!l D74, 

023003 (2006), arXiv:astro-ph/0602538' 

D. Chatterjee and D. Bandyopadhyay, Astrophys. J.| 680, 

686 (2008), arXiv:0712.3171 [astro-ph], 

M. E. Gusakov and E. M. Kant or, ,Phys. Rev.| D78, 



083006 (2008), arXiv:0806.4914 [as tro-phJI 
M. Sinha and D. Bandyopadhyay, [Phys. Rev.| D79, 



123001 (2009), a rXiv:0809.3337 [astro-phj]" 
T. K. Jha, H. Mishra, and V. Sreekanth, Phy s. Revr] C82, 

025803 (2010), arXiv;1002.4253 [hep-ph] 

D. Chatterjee and D. Bandyopadhyay, Phys. Rev.| D75, 

123006 (2007), a rXiv:astro-ph/070 2259 

D. Chatterjee and D. Bandyopadhyay (2007), 

|arXi^:0712.434 7 [astro-ph] [ 

J. Madsen, Ph ys Rev.|D46, 3290 (1992). 

Madsen, Thys. R ev. Lett.] 85, 10 (2000), 

|arXiv :astro-ph/9912418 

M. G . Afford and A. Schmitt, [jTP^s] G34, 67 (2007), 
[arXiv:nucl-th/0608019 



Phys. RevTI 



Phys. Rev. I 



B. A. Sa'd, I. A. Shovkovy, and D. H. Rischke, 1 
D75, 065016 (2007), arXiv:astro-ph/0607643 " 
B. A. Sa'd, 1. A. Shovkovy and D. H. Rischkejl 
D75, 125004 (2007), arXiv:astro-ph/0703016 
M. G. Afford, M. Braby S. Reddy and T. Schafer, 
[Piys. Rev. C75, 055209 (2007), arXiv:nucf-th/0701067 
D. B. Blaschke and J. Berdermann, AIP Conf. Proc. 964. 



290 (2007), arXiv:0710.5243 [hep-pli| 



B. A. Sa'd(2008), |arXiv:0 806.3359 [astro-ph]' 
M. Mannareffi,^ CT Manuef, and B. A. 
Sa'd, Phys. ReO]ett. 101, 241101 (2008), 
,arXiv:0807.3264 [hep-ph], 



[26; 

[27 

[28; 
[29; 
[3o; 

[31 

[32; 

[33; 
[34; 

[35; 

[36 
[37 

[38; 

[39 

[40 

[41 
[42; 

[43; 

[44; 

[45' 
[46; 

[47 

[48; 

[49 

[5o; 



p. Jaikumar, G. Rupak, and A. W. Ste iner, [Phys. Rev.| 

D78, 123007 (2008), arXiv:080 6.1005 [n ucf-th]| 

M. G. Afford, M. Braby , and S. Mahmood gar,|Phys. RevT] 

C81, 025202 (2010), arXiv:0910.2180 [nucf-tt i]| " 

A. Drago, A. Lavagno, and G. Pagfiara, Phys . Rev.| D71, 

103004 (2005), arXiv;astro-ph/0312009 

A. Drago, G. Pagf iara, and I. Parenti, Astrophys. J.| 678, 
L117 (2008), ar Xiv:0704.1510 [astro-ph]] 

L. Lindblom, W. j! Owen, and S. M. Morsink, 
|Phys. Re v. LetEl80, 4843 (1998), a rXiv:gr-qc /9803053. 

B. J. Owen et al, Phys. Rev.| D58, 084020 (1998), 
arXiv:gr-qc/9804044 , 

N. Andersson, U. L Jones, and K. D. Kokko- 
tas, |Mon. Not. Roy Astron. Soc, 337, 1224 (2002), 



L63 (2002), 



Berezhi- 
(2006), 



|arXiv:astro-ph/0111582j 

R. V. Wagoner, Astro physTT] 578, 
|arXiv:a stro-ph/0207589 

A. Drago, G. Pagfiara, and Z. 
ani, Astron. Astrophys. 445, 1053 

arXiv:gr-qc/0405r45r ' 

A. Reisenegger and A. A. Bonacic, Phys . Rev. Lett?] 91, 
201103 (2003), arXiv:astro-ph/0303375 
L. Rezzoffa, F. K. Lamb, and S. L. Shapiro, [Astrophys. J.| 
531, L139 (Mar. 2000). 

L. Rezzoffa, F. K. Lamb, D. Markovic, and S. L. Shapiro, 
Phys. Rev. D64, 104013 (2001). 

L. Rezzoffa, F. K. Lamb, D. Markovic, and S. L. Shapiro, 

Phys. Rev. D64, 104014 (2001)^ 

C. Cuofano and A. Drago, 'J. Phys. Conf. Ser. | 168, 

012008 (2009), arXiv;0903.0349 [astro-ph.HE] , 

C. Cuofano and A. Drago, Phys. Rev. D82, 084027 

(2010), [arXiv:0905. 3368 [astro-ph.HE] 

V. C. A.T%rraro, Astrophys. J. 119, 407 (Mar. 195 4). 

C. Afcock, E. Farhi, and A. Ofinto, [Astrophys. JT] 310, 

261 (1986). 

H. C. Spruit, Astron. Astrophys. 349, 189 (Sep. 1999), 
arXiv:astro-ph/9907138j 

H. C. Spru it, Astron. 'fetrophys.l 381, 923 (Jan. 2002), 
|arXiv:astro- ph/0108207| 

P. Haensel, Nucf. Phys. Proc. Suppf. 24B, 23 (1991). 
J. Braithwaite, Astron. Astrophys. 449, 451 (Apr. 2006), 

arXiv:astro-ph/0509693 

J. Braithwaite and H. C. Spruit, Nature! 431, 819 (Oct. 

2004), arXiv:astro-ph/0502043 

J. Braithwaite and H. C. Spruit, Ast ron. Astrophys. ^ 450, 

1097 (May 2006), |arXiv:astro-ph/051 0287| ' 

J. Braithwaite and^. Nordfund, Astron. Astr ophys. | 450, 

1077 (May 2006), arXiv:astro-ph/0510316 

P. Haensel and J. L. Zdunik, Nuovo Cim.l B121, 1349 

(2006), arXiv:astro-ph/0701258( 



8 



[51] P. Jaikumar , S. R eddy, and A. W. [52] A. Drago and A. Lavagno,' Piys. Lett.| B511, 229 (2001), 
Steiner, Phys. Rev Lett. | 96, 041101 (2006), arXiv:hep-ph/0103209| 
arXiv:nucl-th/0507055| 



